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Abstract 

On 

£NJ ' In this paper we test the hypothesis that the Pioneer anomaly can 

be of gravitational origin by comparing the predicted model-independent 
shifts Aa/o for the semimajor axis of Uranus and Neptune with the 
Voyager 2 radio-technical distance measurements performed at JPL- 
NASA. As in the case of other tests based on different methods and 
data sets (secular perihelion advance, right ascension/declination resid- 
uals over about one century), the orbits of the investigated planets are 
. not affected by any anomalous acceleration like that experienced by 

' the Pioneer 10/11 spacecraft. 
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1 Introduction 

The Pioneer anomaly (Anderson et al. 1998; 2002) consists of an unex- 
pected, almost constant and uniform acceleration directed towards the Sun 

yl Pio = (8.74 ± 1.33) x 10~ 10 m s~ 2 (1) 

detected in the data of both the spacecraft Pioneer 10 (launched in March 
1972) and Pioneer 11 (launched in April 1973) after they passed the thresh- 
old of 20 Astronomical Units (AU; 1 AU is slightly less than the average 
Earth-Sun distance and amounts to about 150 millions kilometers). Latest 
communications with the Pioneer spacecraft, confirming the persistence of 
such an anomalous feature, occurred when they reached 40 AU (Pioneer 11) 
and 70 AU (Pioneer 10). 
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If the Pioneer anomaly is of gravitational origin, it must then fulfil the 
equivalence principle, which is presently tested at a 10~ 12 level (Will 2006) 
and lies at the foundations of the currently accepted theories of gravity. In its 
weak form, it states that different bodies fall with the same accelerations in 
a given external gravitational field. As a consequence, an extra-gravitational 
acceleration like ^4pi should also affect the motion of any other object mov- 
ing, at least, in the region in which the Pioneer anomaly manifested itself. 

In this context, many models have been proposed in order to find some 
possible gravitational explanation of the anomalous acceleration experienced 
by the Pioneer 10/11 spacecraft. E.g., Jaekel and Reynaud (2005) put forth 
a metric linear extension of general relativity which yields an acceleration 
only affecting the radial component of the velocity of a test particle. Brown- 
stein and Moffat (2006) used an explicit, analytical model fitted to all the 
presently available Pioneer 10/11 data points. 

The availability of 

• The latest observational determinations of the secular, i.e. averaged 
over one orbital revolution, extra-advances of perihelia w of the inner 
(Pitjeva 2005a) and of some of the outer (Pitjeva 2006a) planets of the 
Solar System 

• The residuals of the direct observables a cos 5 and 5, where a and 5 are 
the right ascension and the declination, respectively, for the gaseous 
giant planets and Pluto (Pitjeva 2005b) 

processed at the Institute of Applied Astronomy, Russian Academy of Sci- 
ences (IAA, RAS) has recently allowed to 

• Perform clean and unambiguous tests of the possibility that the ac- 
celeration of eq. can also affect the planetary motions in the far 1 
regions of the Solar System (Iorio 2006a; Iorio and Giudice 2006) 

• Dismiss the previously cited mechanisms for the anomalous Pioneer 
behavior (Iorio 2006b; 2006c) 

In this letter we perform a further, independent test of the hypothesis 
that the Pioneer anomaly can be of gravitational origin by exploiting certain 

1 Sanders (2006) disproves the possibility that Api is also present in the inner regions of 
the Solar System by using older data for the perihelion advances of Mercury and Icarus. It 
is straightforward to use the more recent and accurate data for all the inner planets (Pitjeva 
2005a) to put much tighter constraints on any possible anomalous extra-acceleration in 
such a region of the Solar System. 
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short-period, i.e. not averaged over one revolution, features of the semimajor 
axes a of Uranus and Neptune and the radar-ranging distance measurements 
to them performed at Jet Propulsion Laboratory (JPL), NASA, during their 
encounters with the Voyager 2 spacecraft (Anderson et al. 1995). 

The outcome of such a test is consistent with the other ones based on w 
(Iorio 2006a; 2006c) and a cos 5 — 5 (Iorio and Giudice 2006): an acceleration 
like that of eq. (^Q) does not affect the motion of Uranus and Neptune. 

2 The effect of a Pioneer-like acceleration on the 
semimajor axis and comparison with the obser- 
vations 

In (Iorio and Giudice 2006) there are the analytical expressions of the short- 
period shifts induced on the Keplerian orbital elements by a radial, constant 
perturbing acceleration A r , whatever its physical origin may be. For the 
semimajor axis we have 

Aa 2eA r a 2 , ^ , 2e A r . ^ , . , 

— = 7TVl-( cosE - cosE o) = — 7E 9 71~r( cosg -cos£ ), (2) 

where GM is the Kepler's constant, e is the orbital eccentricity, 

< AN >= GM (^) = ;^b (3) 

is the Newtonian acceleration averaged over one orbital period and E is the 
eccentric anomaly which can be expressed in terms of the mean anomaly M 
as (Roy 2005) 

/ e 3\ e 2 3 

E~M+[e sinA^H sin2A4 + -e 3 sin3Al (4) 

\ 8 J 2 8 

The reference epoch is customarily assumed to be J2000, i.e. JD=2451545.0 
in Julian date. From eq. ((21) it can be noted that, whatever the eccentricity 
of the orbit is, 

A " X 0, (5) 



so that Aa/a cannot tell us anything about the impact of an acceleration 
like Api for those planets for which data sets covering at least one full 
orbital revolution exist. To date, only Neptune (P = 164 yr) and Pluto 
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(P = 248 yr) have not yet described a full orbit since modern astronomical 
observations became available after the first decade of 1900. Incidentally, 
let us note that 2 , according to eq. ©, Aa/a = for e = 0. 

The situation is different for Neptune since no secular effects can yet be 
measured for it. Thus, let us use eq. @ and eq. (^Q) for A r getting 



Aa 



= (-2.2882 ± 0.3482) x 10" 6 (cos E - cos E ). (6) 

Nep 



The predicted effect of eq. © can be compared with the latest available 
observational determinations. Pitjeva (2005b) used only optical data (Table 
3 of (Pitjeva 2005b)) for the outer planets (apart from Jupiter) obtaining 
a formal, statistical error 5a = 478532 m for the Neptune's semimajor axis 
(Table 4 of (Pitjeva 2005b)) at JD=2448000.5 epoch (Pitjeva 2006b). By 
re-scaling it by 10 — 30 times in order to get realistic uncertainty we get 



8a 
a 



(optical) 

= (l-3)xl0~ 6 . (7) 

Nep 



It must be compared with eq. © at 3 JD=2448000.5 (E = 107.423 deg) 

(JD = 2448000.5) = (-0.7413 ± 0.1128) x 10~ 6 . (8) 



Aa 

a 



Nep 



Such an effect would be too small to be detected. 

In (Anderson et al. 1995) the radio-technical data of the Voyager 2 
encounter with Neptune were used yielding a unique ranging measurement 
of a (Julian Date JD=2447763.67); eq. ©, evaluated at such epoch {E = 
106.012 deg), predicts 



Aa 

a 



(JD = 2447763.67) = (-0.7954 ± 0.1210) x 10~ 6 . (9) 

Nep 



By assuming for Aa the residuals with respect to the DE200 JPL ephemerides 
used in Table 1 of (Anderson et al. 1995), i.e. 8224.0 ± 1 km, one gets 



Aa 

a 



(ranging) 

= (1.8282 ± 0.0002) x 10~ 6 . (10) 

Nep 



2 In the circular orbit limit, Anderson et al. (1998; 2002) use the erroneous formula 
Aa/a — — At/An and apply it to Mars and the Earth to show that an extra-acceleration 
like Api cannot exist in that regions of the Solar System. 

3 For Neptune E = 128.571 deg at JD=2451545.0. 
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This clearly rules out the prediction of eq. ©. 

The same analysis can also be repeated for Uranus (P = 84.07 yr) for 
which no modern data covering a full orbital revolution were available at the 
time of the Anderson et al. (1995) work; as for Neptune, one radar-ranging 
distance measurement is available from the Voyager 2 flyby with Uranus 
(JD=2446455.25). The prediction of eq. ©, with eq. (UJ for A r , for the 
flyby epoch 4 (E = 8.860 deg) is 



Aa 

a 



(JD = 2446455.25) = (-3.3576 ± 0.5109) x 10 -6 . (11) 

Ura 



Table 1 of (Anderson et al. 1995) yields for the DE200 residuals of the 
Uranus' semimajor axis Aa = 147.3 ± 1 km, so that 



Aa 

a 



(ranging) 

(JD = 2446455.25) = (0.0513 ± 0.0003) x 10~ 6 . (12) 

Ura 



Also in this case, the effect which would be induced by Api Q on Aa/a is 
absent. 

It maybe interesting to note that the paper by Anderson et al. (1995) 
has been used as a basis for other tests with the outer planets using different 
methods. E.g., Wright (2003) and Sanders (2006) adopt the third Kepler's 
law. Basically, the line of reasoning is as follows. In the circular orbit limit, 
let us write, in general, P = 2ira/v; in particular, the third Kepler law 
states that P = lix^J a 3 /K p , where K p = GMq. If we assume that K p may 
vary by AK p for some reasons 5 inducing a change in the orbital speed, then 
Av/v = (1/2) AK p / K p . In general, for an additional radial acceleration 
acting upon a test particle in circular orbit AA, AA/A = 2 Av/v: thus, we 
have 

Now, a measurement of the planet's velocity is needed to get AK p /K p (or, 
equivalently, AA/A): since v = na, where n is the orbital frequency, this 
requires a measurement of both a and n, while in our case we only use a. 
Moreover, the measurement of the orbital frequency pose problems for such 
planets which have not yet completed a full orbital revolution, as it was the 
case for Uranus and Neptune at the time of the analysis by Anderson et al. 
(1995). For Neptune, according to the last row of Table 2 of (Anderson et 

4 For Uranus E = 70.587 deg at JD=2451545.0. 
5 E.g. due to dark matter (Anderson et al. 1995). 
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al. 1995), AK P /K™ eas = (-2.0 ± 1.8) x 10~ 6 , while A Pio /A N = (-133.2 ± 
20.3) x 10~ 6 . As can be noted, also in this case, the answer is negative but 
the accuracy is far worse than in our test. 

3 Discussion and conclusions 

In this paper we have used the NASA-JPL radio-technical ranging mea- 
surements to Uranus and Neptune performed during the Voyager 2 flybies 
(Anderson et al. 1995) in order to make a model-independent test of the 
hypothesis that an anomalous acceleration of gravitational origin like that 
detected in the data of the Pioneer 10/11 spacecraft may also affect the or- 
bital motion of such planets. The answer is neatly negative, as in previous 
tests involving the perihelion secular advance of Uranus (Iorio 2006a; 2006c) 
and the right ascension/declination residuals of Uranus, Neptune and Pluto 
over about one century (Iorio and Giudice 2006). 

Thus, in regard to the celestial bodies lying at the edge of the region 
in which the Pioneer anomaly manifested itself (~ 20 — 70 AU), or entirely 
residing in it, the present-day situation can be summarized as follows 

• Uranus (a = 19.19 AU). 3 model-independent tests 

— Secular advance of perihelion (almost one century of optical data 
processed at IAA, RAS): negative 

— Right ascension/declination residuals (almost one century of op- 
tical data processed at IAA, RAS): negative 

— Short-period semimajor axis shift (1 radar-ranging measurement 
at epoch JD=2446455.25 by JPL, NASA): negative 

• Neptune (a = 30.08 AU). 2 independent tests 

— Right ascension/declination residuals (almost one century of op- 
tical data processed at IAA, RAS): negative 

— Short-period semimajor axis shift (1 radar-ranging measurement 
at epoch JD=2447763.67 by JPL, NASA): negative 

• Pluto (a = 39.48 AU). 1 test 

— Right ascension/declination residuals (almost one century of op- 
tical data processed at IAA, RAS): negative 
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In all such tests the observationally determined quantities— obtained at JPL 
and IAA independently and without having the Pioneer anomaly in mind at 
all— have been compared to unambiguous theoretical predictions based on 
the effects induced by a radial, constant and uniform acceleration with the 
same magnitude of that experienced by Pioneer 10/11, without making any 
assumptions about its physical origin. 

In addition, we may also consider the perihelion-based negative tests for 
Jupiter (a = 5.20 AU) and Saturn (a = 9.53 AU) (Iorio 2006c), based on the 
model by Brownstein and Moffat (2006) fitted to all the presently available 
data points of Pioneer 10/11. 

In conclusion, it seems more and more difficult to accept the possibility 
that some modifications of the current laws of Newton-Einstein gravity may 
be the cause of the Pioneer anomaly. 
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